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ABSTRACT 


The  eain  in  sipnal-to-noise  ratio  with  resnect  to  a bandnass 
filter  was  investittated  as  a function  of  imut  siunai-to -noise  ratio  for  a Wiener 
filter,  a matched  filter,  and  the  three  comnonent  adaDti\e  filter.  Each  pro- 
cessor  was  annlied  to  two  southern  Kamchatka  events  recorded  at  a single 
site  at  the  Alaskan  Long  Period  Array  (ALPA)  and  buried  in  two  noise  sam- 
nles  recorded  there. 

A r enr esentative  southern  Kamchatka  earthquake  seismogram 
was  formed  by  averaging  snectra  of  events  from  that  region.  The  Wiener 
filter  used  the  resulting  average  oower  snectrum  but  no  signal-noise  correla- 
tion terms  to  achieve  a i to  6 dB  maximum  gain  over  the  bandnass  filter.  The 
average  event  was  also  used  as  a matched  filter,  which  had  a gain  ol  mom  6 
to  10  dB  over  the  bandnass  filter.  With  the  heln  of  a Monte  Carlo  simulation 
to  determine  the  nest  form  for  the  filter  weights,  the  three  comnonent  adan- 
tive  filter  had  a gain  of  from  10  to  17  dB  over  the  bandnass  filter. 
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SECTION  I 


INTRODUCTION 


This  renort  describes  the  performance  of  three  important  pro- 
cessors for  seismic  data,  the  Wiener  filter,  the  matched  filter,  and  the  three 
component  adaptive  (TCA)  processor.  Each  has  been  optimized  as  far  as  pos- 
sible for  the  present  application. 

The  data  used  here  were  single-site  Alaskan  Long  Period  Array 
(ALPA)  signals  and  noise  samples.  Single-site  data  were  used  because  the  de- 
tection problem  is  more  severe  for  a single  instrument  than  for  an  array,  and 
because  the  trend  at  present  is  toward  deployment  of  single-site  instrumenta- 
tion (Seismic  Research  Observatories)  and  small  arrays  (Iranian  Long  Period 
Array).  ALPA  data  were  used  because  they  are  of  high  quality  and  because 
the  nature  of  noise  and  signals  at  ALPA  is  well  known. 

Synthetic  data  were  prepared  by  selecting  a Z048-point  segment 
of  seismic  noise,  sampled  every  two  seconds,  and  adding  to  it  a suitably  scaled 
signal,  originally  recorded  at  high  signal-to-noise  ratio  (SNR),  to  form  a com- 
posite trace.  The  ratio  of  the  scaled  signal  peak  to  the  root-mean-square 
(RMS)  value  of  the  noise  preceding  the  time  at  which  the  signal  was  added  was 
the  true  input  signal-to-noise  ratio.  For  the  matched  and  Wiener  filters,  dif- 
ferent scale  factors  were  used  for  each  trace  so  that  all  had  the  same  SNR. 

For  the  TCA,  a common  scale  factor  was  used  so  that  the  signal  was  not  dis- 
torted. This  procedure  has  the  advantage  over  use  of  signals  as  recorded  that 
the  true  SNR  and  peak  arrival  time  are  known  exactly,  and  can  be  varied  by 


the  experimenter 


The  signal  buried  in  noise  is  called  the  test  event.  Another 
event  is  required  for  implementation  of  the  Wiener  and  matched  filter  and  is 


cnlled  the  reference  event.  Descriptions  of  the  noise  samples  and  test  events 
used  in  this  study  are  given  in  Table  1-1,  and  the  events  are  plotted  in  Figure 
1-1. 

The  general  experimental  procedure  was  as  follows  The  test 
event  was  buried  in  noise  at  various  known  SNR's,  and  the  nrocessor  under 
study  applied  to  the  composite  trace,  The  ratio  of  the  largest  peak  in  the  sig- 
nal arrival  interval  to  the  RMS  value  of  the  preceding  noise  was  calculated. 
The  SNR  was  also  calculated  for  the  same  composite  trace  filtered  between 
0.024  Hz  and  0.  O'iO  Hz,  the  frequency  range  where  the  SNR  of  seismic  events 
is  expected  to  be  large.  The  difference  (in  dB)  between  these  ratios  was  the 
gain  of  the  processor  over  the  bandpass  filter  and  was  plotted  as  a function  of 
input  SNR. 

It  is  useful  to  present  gains  rather  than  absolute  SNR's  because 
gains  do  not  increase  indefinitely  with  input  SNR,  but  saturate  at  some  maxi- 
mum value.  The  bandpass  filter  was  chosen,  as  a standard  because  it  is  a 
routine  form  of  signal  improvement  whose  performance  is  well  known. 

A schematic  form  for  the  results  is  shown  in  Figure  1-2.  The 
horizontal  axis  is  input  SNR  and  the  vertical  axis  is  SNR  gain  over  the  band- 
pass filter.  Both  scales  are  in  dB,  although  no  values  are  shown. 

At  low  values  of  input  SNR  the  signal  is  so  far  below  the  noise 
that  the  peak  of  the  processor  output  is  not  contributed  by  the  signal.  Conse- 
quently, the  gain  is  not  a function  of  the  input  SNR.  At  sufficient  input  SNR, 
the  nrocessor  output  is  determined  by  the  signal,  although,  hopefully,  the 
bandpass  filters'  output  is  still  due  to  noise.  At  this  point  the  gain  of  the 
processor  begins  to  increase.  At  some  higher  SNR  the  bandpass  filter  peak 
is  due  to  the  signal,  and  thereafter  little  gain  is  oiitained.  In  some  cases  the 
gain  may  decrease.  In  any  case,  the  SNR  gain  reaches  a maximum,  as  shown 
in  Figure  1-2,  and  this  maximum  gain  is  taken  here  as  a measure  of  nroces- 
sor performance. 
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FIGURE  1-2 

SCHEMATIC  GAIN  VERSUS  INPUT  SNR 
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ders,  or  occasionally  no  shoulders  at  all.  Since  the  exact  behavior  of  both 
signal  and  noise  are  known  in  a study  such  as  this,  each  unusual  feature  could 
be  exolained  by  a detailed  examination  of  the  data  and  the  processor.  Such  an 
exnianation  is  not  attempted  because,  while  the  basic  process  giving  rise  to 
signal  and  noise  are  deterministic,  from  the  geophysicist's  point  of  view  they 
are  random.  Both  noise  and  signal  are  thought  of  as  being  drawn  from  some 
poDulation,  and  are  unlikely  to  occur  again.  Thus  only  the  most  general  fea- 
tures of  the  results  are  repeatable,  and  deviations  from  them  must  be  expect- 
ed and  treated  as  fluctuations. 
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SECTIOM  11 


WIENER  FILTER 


I 

I 


The  Wiener  filter  is  one  whose  output  is  closest  to  the  desired 
signal  in  the  least  mean  square  sense.  In  the  frequency  doniain,  it  takes  the 
following  form: 


F = 


(1>  + <l> 

s s s n 


(/)  + (I)  + (/)  + <l> 

ss  sn  ns  nn 


(Il-l) 


where  the  d>'s  are  the  cross-  and  auto-power  spectra  of  the  signal  and  noise. 
Thus,  if  X is  the  test  event  which  it  is  desired  to  extract  from  the  noise  sam- 
ple z , the  <b's  are  defined  as  follows; 


(1)  = 
na 

z(o>) 

(h  ~ 

sn 

x{l,)) 

• z':i:(co) 

(11-3) 

(/)  = 
ns 

z((o) 

• x*(w) 

(11-4) 

(l> 

x(m) 

• X^  ( M ) 

(11-5) 

where  o>  is  the  angular  frequency,  the  symbol  ~ denotes  a Fourier  trans- 
form, and  * denotes  a complex  conjugate. 

The  Wiener  filter  in  the  form  of  equation  (71-1)  is  in  a sense 
trivial.  Rewriting  it  using  equations  (ll-Z)  through  (11-S),  we  have  the  follow- 
ing ; 


F = 


~(o)  )x-!=(  (o)  -1  x(dJ  )7*(  (rl)  1 7'(  (O  )X*(  e)  ) + 7.  (u>)7.^(ui) 


(IJ-b) 
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The  denominator  of  equation  (II-6)  can  be  factored  as  can  the  numerator,  as 
follows : 


!■ 


1 


F = 


x((,j)  (x((,jJ  + z(m)) 
(x(oj)  + z(w))  (x((o)  + 'z(oj))''' 


(11-71 


The  output  of  the  filter  is  F times  the  input  x(ai)  + z'(oj)  , and  is: 


|x(w)  +?((o)| 


x((ti)  (x((a)  +z{(tj ) ) 


(x((u)  + z:(oj))  (x(w)  + z(oj)) 


|x(oj)+~(c4l  =x(<o).  (11-8) 


I'his  result  means  that  if  the  signal  and  noise  are  known  exactly,  the  signal 
can  be  extracted  exactly  from  the  noise.  The  utility  of  the  Wiener  filter  is, 
of  course,  that  the  form  of  equation  (11-1)  can  also  be  used  when  signal  and 
noise  are  only  approximately  known. 

The  problem,  then,  is  to  estimate  the  terms  in  equation  (11-1), 
The  noise  power  spectrum  can  be  found  from  the  noise  preceding  the 

signal  arrival,  and,  as  usual,  we  will  assume  that  the  noise  is  stationary  so 
that  this  power  spectrum  can  be  used  to  estimate  the  power  spectrum  in  the 
presence  of  the  signal.  A model  for  the  signal  must  be  assumed  in  order  to 
calculate  '('gg«  A model  is  also  required  to  find  the  crosspower  terms 
and  no  satisfactory  way  of  doing  this  has  been  found. 

The  noise  power  spectrum  was  estimated  by  segmenting  a 
2048-point  series  into  8 sections  of  256  points  each,  filling  with  zeros  to  make 
2048  points,  and  Fourier  transforming.  At  each  frequency,  the  power  was  the 
sum  of  the  squares  of  the  real  and  imaginary  parts.  These  powers  were  aver- 
aced  at  each  frequency  over  all  8 segments  to  increase  the  reliability  of  the 
estimates. 

The  signal  power  spectrum  was  fovind  by  averaging  power  spec- 
tra from  seven  southern  Kamchatka  events  together.  Lach  had  been  normalized 


11-2 


to  unit  area  before  averauinu,  to  eliminate  the  influence  of  magnitude 
average  power  spectra  are  shown  in  Figure  11-1. 


The 


j Phase  spectra  for  these  events  were  corrected  for  jumps  of 

12—  by  adding  or  subtracting  2.  rr  whenever  such  a jump  was  found.  Then 
their  linear  trend  was  removed  by  subtracting  the  quantity 
i phase  angle  at  frequency  f , where  f^  is  the  Nyquist  frequency  and  was 

I the  phase  at  the  Myquist  frequency  of  the  transverse  component.  This  had  the 

effect  of  moving  each  signal  without  distortion  so  that  it  was  centered  on  z.ero 
time. 


Finally,  the  seven  corrected  phase  spectra  were  averaged  to- 
gether for  each  component  and  are  shown  in  Figure  II-Z.  In  regions  where 
the  signal  power  is  low,  the  phase  spectrum  is  meaningless,  but  over  the  sig- 
nal band  the  slope  of  the  phase  spectrum  is  indicative  of  the  phase  velocity  of 
the  signal. 

An  amplitude  spectrum  was  calculated  for  each  component  by 
extracting  the  square  root  of  the  power  spectrum,  and  the  inverse  transform 
of  this  amplitude  and  phase  spectrum  found  and  plotted  in  Figure  11-3.  These 
synthetic  signals  show  the  effects  of  multipathing  (which  is  also  visible  in  the 
scalloping  of  the  power  spectra  in  Figure  II- 1),  of  dispersion,  and  of  the  higher 
group  velocity  of  Fove  waves  than  Rayleigh  waves.  This  time-domain  signal 
will  also  be  used  as  a matched  filter,  and  is  denoted  as  s(t). 

The  average  correlation  coefficient  between  the  synthetic  sig- 
nal found  here  and  seven  independent  events  from  southern  Kamchatka  was 
0.69,  with  a standard  deviation  of  0.08.  The  average  correlation  between 
KAM-2  31-19AL  and  the  same  seven  events  was  0.64,  with  a standard  de\-ia- 
tion  of  0.  18.  These  mean  correlations  are  indistinguishable,  but  the  stand- 
ard deviation  for  the  synthetic  signal  is  less  thati  for  the  ri'al  one.  This 
means  that  although  there  is  no  advantage  in  terms  of  gain  in  using  the 


i 

i 

1 

i 


lU  •> 


T f.T  r s <■ 


H .1  fii.ll 


I 

0.  i5 


F1GUFU-.  II -1 

powi;r  s''>[;ctra  for  sYNTniiiic  kvknt 


synthetic  waveform,  there  is  a hiiiher  likelihood  of  detectin;^  a sij^nal  usinj4 
the  synthetic  waveform,  since  detection  is  a yes-no  decision  depending  on  a 
th  reshold 

The  assumption  is  ordinarily  made  that  signal  and  noise  are 
uncorrelated  so  that  terms  like  may  be  ignored.  VVhile  this  is  true  on 

the  average,  over  the  short  sections  of  data  used  here  there  will  be  some  ap- 
parent correlation  between  signal  and  noise  due  to  the  random  nature  of  the 
noise.  Furthermore,  at  low  enough  values  of  signal- to- noise  ratio,  the  cor- 
relation term  (/>  will  dominate  the  numerator  of  eciuation  (11-1),  since  it  is 
s n 

proportional  to  the  first  power  of  the  signal  amplitude,  while  </)^^  is  propor- 
tional to  the  square  of  the  signal  amplitude.,  Consequently,  it  is  important  to 

include  (t>  in  the  Wiener  filter, 
s n 

The  terms  contributing  to  the  Wiener  filter  for  the  transverse 
component  are  plotted  in  Figure  11-4.  Here  a signal-to - nois e ratio  of  1 0 dB 
has  been  used,  and  <l>  , ,f>  , and  </>  plotted  on  the  same  scale  as  a func- 

ss'  • nn  I 

tion  of  frequency.  Since  the  phase  angle  of  the  noise  is  random,  the  phase  of 

(f)  is  also  random,  and  is  not  shown, 
sn 

The  influence  of  the  signal-noise  correlation  terms  is  shown 
explicitly  in  Figure  11-^.  where  the  Wiener  filter  is  shown  for  the  same  sig- 
nal and  noise,  at  10  dB  s ignal-to -noise  ratio,  with  and  without  the  correlation 
terms.  For  comparison,  the  bandpass  filter  used  here  is  also  shown.  It  is 
clear  that  their  inclusion  makes  a significant  difference  to  the  filter  weights. 

The  term  </»  in  the  numerator  also  introduces  a phase  shift  into  the  filter, 
sn 

Having  established  the  importance  of  the  signal-noise  correla- 
tion terms,  the  next  step  is  to  calculate  them  for  an  unknown  signal.  The 

average  amplitude  of  </»  is  easy  enough  to  find,  but  its  phase  is  unpredict - 

sn 

able,  since  the  noise  phase  is  a random  variable  depending  on  the  details  ot 
the  time  domain  motion,  and  can  take  on  with  equal  prol)al)ility  any  value  over 
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FIGURE  II-4 

AND|0j.r^l  FOR  TRANSVERSE  COMPONENT 
OF  REFERENCE  EVENT  AT  SNR  = 10  dB 


FIGURE  II-5 

WIENER  FILTER  WITHOUT  (a)  AND  WITH 
(b)  SIGNAL- NOISE  CORRELATION 
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its  whole  range.  The  probability  distribution  function  of  the  power,  by  con- 
trast, is  peaked  at  the  mean  power,  so  that  the  power  at  a given  frequency  is 
much  more  predictable  than  the  phase.  There  are  many  time  series  having 
the  same  power  spectrum  but  different  phase  spectra,  and  these  time  series 
differ  in  that  their  energy  appears  at  different  times.  In  practice,  this  ran- 
domness in  the  distribution  of  energy  in  the  time  domain,  and  the  resulting 
randomness  of  the  phase  spectrum,  was  impossible  to  model. 

To  examine  the  performance  of  the  Wiener  filter,  it  was  cal- 
culated for  a number  of  signal-to-noise  ratios  and  applied  to  signals  buried 
at  the  appropriate  level  in  noise.  Various  models  of  the  noise  phase  were 
examined.  The  noise  phase  was  taken  to  be  its  expected  value  of  zero,  or  to 
be  a random  number,  or  to  be  equal  to  the  signal  phase.  None  of  these  choices 
was  as  effective  as  omitting  the  correlation  term  entirely.  Consequently,  all 
results  for  the  Wiener  filter  presented  here  assume  that  <t>^^  = 0.  The  filter 
weights  are  then  as  follows: 


F 


(h  + (/) 
s s nn 


(II-9) 


The  gain  of  the  Wiener  filter  over  the  bandpass  filter  is  shown 
as  a function  of  input  SNR  in  Figures  11-6  and  11-7.  In  the  first  figure,  the 
noise  sample  is  varied  while  the  test  event  is  held  constant,  and  the  opposite 
is  true  in  Figure  II-7.  More  correlation  is  seen  between  results  with  a com- 
mon test  event  than  those  with  a common  noise  sample,  because  the  perfor- 
mance of  ihe  filter  is  principally  affected  by  the  correlation  between  test  and 
reference  events  rather  than  by  the  uniformly  low  correlation  between  the 
noise  and  reference  event. 

Gains  are  rather  low  compared  to  those  from  other  processors 
studied  here.  Peak  gains  varied  from  2 to  6 dR  and  occurred  at  from  6 to  16 
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FIGURE  11-6 

GAIN  VERSUS  INPUT  SNR  FOR 
WIENER  FILTER  WITH  CONSTANT  SIGNAL 
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Gain  Over  Bandpass  Filter  (dB) 


Input  SNR  (dB) 


FIGURE  II-7 

GAIN  VERSUS  INPUT  SNR  FOR 
WIENER  FILTER  WITH  GONSTANT  NOISE 


11-12 


ik  6 


r 

‘.1 
U! 

5 dB  input  SNR.  Output  waveforms  for  an  input  SNR  of  1 0 dB  are  shown  for  the 

I:  Wiener  filter  in  Figure  11-8  and  for  the  bandpass  filter  in  Figure  11-9,  for 

S'  K AM-2 3 1 - 1 9 AL,  buried  in  noise  sample  NCI- 3 1 1 - 08 AL, . 

II 

[Without  the  signal-noise  correlation  terms  the  Wiener  filter  is 

just  a zero-phase  filter  with  a shape  which  suppresses  more  of  the  noise  (and 
of  the  signal)  than  does  the  simple  rectangular  bandpass  filter.  Since  no  qual- 
itatively different  operations  are  performed,  it  is  not  surprising  that  the  gain 
over  the  bandpass  filter  should  not  be  large.  We  may  also  conclude  from  the 
small  changes  in  performance  when  the  noise  sample  was  changed  that  the 
noise  and  signal  spectral  peaks  are  fairly  well  separated,  at  least  for  that  sig- 
nal. 
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SECTION  III 

MATCHED  FILTER 


Prewhitening  followed  by  matched  filtering  is  the  optimum  de- 
tection scheme  (Capon,  Greenfield,  and  Lacoss,  19£>9).  Prewhitening  is  ac- 
complished by  dividing  the  transform  of  the  trace  to  be  processed  by  the  pow- 
er snectrum  p((o)  of  the  noise.  The  matched  filter  used  here  is  the  synthet- 
ic signal  dev'eloped  in  the  previous  section  on  Wiener  filtering,  and  is  con- 
volved with  the  prewhitened  input  trace  to  yield  an  output  y . 

Using  the  notation  of  the  previous  section,  where  x(t)  is  the 
test  event,  z(t)  the  noise,  and  s(t)  the  matched  filter  waveform,  the  opera- 
tions of  matched  filtering  and  prewhitening  can  be  combined  into  one  expres- 
sion in  the  frequency  domain,  as  follows: 


x((ij)  + z(a»)  s=''={a») 
p(w) 


(III-l) 


Comparing  this  expression  to  that  for  the  Wiener  filter  without 
crosscorrelation  terms  (equation  (II-9)),  it  is  evident  that  the  Wiener  filter 
outnut  is  approximately  the  prewhitened  matched  filter  output  convolved  with 
the  reference  waveform.  Intuitively,  it  would  seem  that  such  convolution  can 
only  reduce  the  signal-to-noise  ratio  and  that  the  Wiener  filter  should  display- 
less  gain  than  does  the  matched  filter.  It  will  be  seen  that  this  is  the  case. 

The  output  of  the  filter  (equation  (111-1))  is  proportional  to  the 
input,  but  with  unknown  scale  factor.  In  order  to  find  correct  peak  amplitudes, 
we  calculate  a multiplicative  scale  factor  for  the  reference  waveform  that  will 
ensure  that  matched  filtering  a waveform  with  itself  in  the  absence  of  noise 


lll-l 


will  yield  a peak  amplitude  equal  to  that  of  the  original  signal.  This  require- 
ment is  satisfied  if: 


g s[m)  s*((0) 


p(w) 


= s(t) 


max 


(111-2) 


where  g is  the  desired  scale  factor,  1 denotes  inverse  Fourier  transforma- 
tion, and  the  subscript  'max'  denotes  taking  the  maximum  of  a waveform. 
Then: 


s(t) 


g = 


(111-3) 


s'(g))  s'* (to) 
p(o>) 


max 


The  function  in  the  denominator  of  equation  (111-3)  is  the  prewhitened  autocor- 
relation function  of  the  reference  waveform. 

The  scale  factor  g thus  found  is  not  a function  of  the  test 
event,  so  that  it  can  be  used  to  extract  an  unknown  event.  The  matched  filter 
peak  will  be  somewhat  smaller  than  the  true  value  because  of  less  than  per- 
fect correlation  between  test  and  reference  events,  but  this  loss  will  be  only 
about  3 dB  for  a correlation  of  0.7, 

Gains  as  a function  of  input  SNR  are  shown  in  Figures  Ill-l 
and  II1-2  in  the  same  format  as  in  the  previous  chanter.  Figure  IIl-l,  for 
constant  signal,  shows  some  similarity  in  the  shape  of  the  gain  curves  as  the 
noise  is  varied,  as  for  the  Wiener  filter,  although  the  levels  are  different. 
Again,  this  is  because  the  gain  is  determined  by  the  test  event-reference 
event  correlation,  reference  event-noise  correlation  being  uniformly  low. 

There  is  no  obvious  similarity  between  gains  for  tlie  same 
component  from  different  e\'ents  in  Figure  III-2,  where  the  noise  sample  is 
held  constant.  This  is  because  the  most  important  factor,  the  correlation 
between  test  and  reference  event,  is  not  constant.  Furthermore,  prewbiten- 
ing  removes  some  of  the  effect  of  \ arying  the  noise  sample. 
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FIG  UR  f:  lll-l 

GAIN  VF'.RSUS  INPUT  SNR  FOR 
MATCHFD  FILTFR  WITH  CONSTANT  SIGNAL 
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There  is  clear  similarity  between  Figures  IIl-l  and  1II-2  with 


Figures  11-6  and  11-7,  the  corresponding  Wiener  filter  results.  The  Wiener 
filter  results  are  roughly  parallel  to  the  matched  filter  results,  component 
by  component,  but  lie  about  3 dB  below  them.  This  is  consistent  with  the  dis- 
cussion earlier  in  this  chapter  where  it  was  pointed  out  that  the  Wiener  filter 
output  is  just  the  matched  filter  output  convolved  with  a reference  waveform. 
This  process  of  convolution  can  now  be  seen  to  lose  about  3 dB  in  terms  of 
output  SNR. 

A plot  of  matched  filter  outputs  for  K AM-2 3 1 - 1 9 ALi  buried  at  a 
true  SNR  of  10  dB  in  noise  sample  NOl-31  1 -OSAL  is  shown  in  Figure  111-3. 

The  peaks  occur  at  the  signal  arrival  time  rather  than  at  the  peak  arrival  time 
as  do  processors  which  do  not  concentrate  the  signal  energy  into  a single  point. 
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SECTION  IV 

THREE  COMPONENT  ADAPTIVE  FILTER 


The  three  component  adaptive  filter  (TCA)  examined  here  has 
been  described  by  Strauss  (1976).  It  makes  use  of  all  three  components  of 
the  seismic  record  simultaneously  and  designs  a data-dependent  filter  for 
each  of  a large  number  of  overlapping  segments  of  the  record.  The  filter  ap- 
plied to  the  vertical  and  radial  components  has  large  filter  weights  at  frequen- 
cies where  the  phase  angle  between  those  components  is  90  degrees;  i.  e.  , 
Rayleigh  waves  only  are  passed.  The  transverse  filter  coefficients  are  large 
when  transversely-polarized  energy  is  present,  so  that  Eove  waves  are  pass- 
ed. Only  energy  in  phase  with  the  vertical  motion  is  considered  in  the  design 
of  the  transverse  filter,  in  order  to  avoid  confus ion  of  Rayleigh  waves  with  off- 
azimuth  Love  waves. 

The  basic  design  of  this  filter  is  fixed,  and  optimization  of  its 
performance  can  come  only  through  variation  of  parameters.  The  parameters 
examined  here  were  the  length  of  the  overlapping  data  segments,  the  number 
of  segments  overlapped  at  each  time  point,  and  the  dependence  of  the  filter 
weights  on  the  data. 

It  was  found  that  there  was  no  significant  difference  in  perfor- 
mance when  64  point  or  128  point  overlapping  segments  were  used,  and  this 
result  is  consistent  with  previous  studies.  Some  degradation  was  observed 
when  32  point  and  256  point  segments  were  used,  so  that  a 64  point  segment 
was  used  in  this  study  in  the  interest  of  reduced  computational  tinie. 

Next,  four  segments  of  data  were  overlapped,  rather  than  only 
two  as  in  all  previous  studies.  It  was  hoped  that  this  would  result  in  more 


i 
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rapid  adaptation,  due  to  shorter  intervals  between  initial  segment  points,  but 
without  reduction  in  frequency  resolution  due  to  decreased  segment  length. 
However,  no  difference  in  performance  between  this  processor  and  that  over- 
lapping only  two  segments  was  found,  although  from  3Z  to  Z56  point  segments 
were  examined.  Therefore,  the  original  version  was  retained,  again  in  the 
interest  of  computational  speed 


In  both  Love  and  Rayleigh  filters,  the  filter  weight  is  a function 
of  an  angle.  For  the  Rayleigh  filter  the  relevant  angle  is  that  between  the 
phases  of  the  vertical  and  radial  Fourier  components  of  motion,  while  for  the 
Love  filter  it  is  the  angle  in  space  between  the  horizontal  particle  motion  and 
the  direction  of  propagation.  Since  it  would  be  impossible  to  examine  all  func- 
tional foi’-ns  for  these  filter  weights,  a Monte  Carlo  method  was  adopted  to 
find  reasona.,-vi  candidates. 


The  model  used  for  the  Love  wave  filter  design  is  illustrated 
in  Figure  IV-1.  An  on-azimuth  Love  wave  and  noise  are  present.  The 
noise  is  assumed  to  have  equal  amplitude  N in  transverse  and  radial  direc- 
tions and  random  phase  with  respect  to  the  vertical  Fourier  component  of 
motion.  Since  only  motion  in  phase  with  the  vertical  is  considered,  the  ef- 
fective noise  amplitude  is  its  true  amplitude  multiplied  by  the  cosine  of  its 
(random)  phase  angle  7 , The  angle  j3  in  Figure  IV-1  is; 


j8  = tan 


-1 


N cos  y 


S + N cos  y 
t t 
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and  we  desire  the  distribution  of  this  angle,  given  that  the  phase  angles  of  the 
transverse  and  radial  components  of  the  noise,  7^  and  7 , are  distributed 
uniformly  over  the  range  of  0 to  Z“. 

This  problem  would  be  difficult  to  solve  analy tically,  but  lends 
itself  to  the  Monte  Carlo  method.  At  each  signal-to-noise  ratio  of  interest, 
we  choose  a large  number  (say  1000)  of  pairs  of  independent  numbers  7 


FIGURE  lV-1 

LOVE  WAVE  FILTER  DESIGN 
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and  y uniformly  distributed  between  0 and  2~,  and  compute  the  correspond- 
r 

inc  values  of  /3  . The  probability  distribution  function  is  found  by  counting  the 
number  of  times  each  value  of  /3  occurs. 

The  signal-to -noise  ratio  used  here  is  a narrowband  ratio  of 
two  Fourier  amplitudes,  rather  than  the  wideband  ratio  of  peak  to  root-mean- 
square  amplitude  used  to  measure  gain.  The  narrowband  ratio  is  strongly 
dependent  on  the  distribution  of  signal  and  noise  energy  with  frequency.  The 
phase  of  the  Love  wave  is  arbitrary  with  resoect  to  that  of  the  vertical  Ray- 
leigh so  the  effective  Love  wave  amplitude  is  smaller  than  that  found  from  a 
spectrum. 

We  can  get  an  idea  of  the  best  weighting  function  by  calculating 
the  distribution  function  of  /S  for  a few  values  of  the  signal-to-noise  ratio  in 
the  range  of  interest.  Such  functions  are  shown  in  Figure  IV-2  for  narrow- 
band  signal-to-noise  ratios  of  4,  2,  1,  0.5,  and  0,25,  corresponding  to  wide- 
band signal-to-noise  ratios  from  about  20  to  -4  dB  for  a typical  test  event  and 
noise  sample  used  here. 

The  weighting  function  was  chosen  to  be  the  one  shown  in  Fig- 
ure IV-3  with  center  point  d = 0.  0°  for  frequencies  between  0.  024  Hz  and 

o 

0.059  Hz  and  zero  outside  that  passband.  The  probability  of  getting  any  use- 
ful information  at  frequencies  where  the  signal-to-noise  ratio  is  less  than 
about  one  can  be  seen  to  be  very  small  from  Figure  IV-2  so  these  frequencies 
outside  the  signal  band  are  ignored. 

The  Rayleigh  wave  filter  was  designed  in  much  the  same  way. 
Figure  IV-4  shows  the  relevant  angles,  which  are  now  in  phase  space  rather 
than  coordinate  space  as  for  the  Love  wave  filter.  The  equation  for  the  angle 
>1'  , which  is  just  90  degrees  for  pure  Rayleigh  waves,  is 
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and  a Monte  Carlo  simulation  gave  the  distribution  functions  shown  in  Figure 
IV-5.  Again,  Fourier  comoonents  outside  the  usual  passband  were  zeroed, 
and  inside  the  passband  the  weighting  function  shown  in  Figure  IV-3  was  ap- 
plied with  90°. 

Three  combinations  of  signal  and  noise  samples  were  process- 
ed, as  usual.  The  results  are  nresented  in  Figures  lV-6  for  constant  signal 
and  Figure  lV-7  for  constant  noise.  As  for  the  Wiener  and  matched  filters, 
there  is  stronger  correlation  between  curves  for  constant  signal  than  for  ones 
with  constant  noise.  Presumably  this  is  because  one  signal  fits  the  model 
used  to  find  the  filter  weights  better  than  another,  whereas  all  noise  samples 
agree  about  the  same.  Output  waveforms  for  the  TCA  at  input  SNR  of  1 0 dB 
are  shown  in  Figure  IV-8.  The  effect  of  segmenting  the  data  is  clear. 

Gains  for  the  TCA  varied  from  9 to  17  dB,  and  occurred  over 
a range  of  input  SNR  from  8 to  16  dB.  These  gains  are  greater  than  those  ob- 
tained from  the  matched  or  Wiener  filters.  It  may  be  that  this  superior  per- 
formance is  due  to  the  fact  that  the  TCA  needs  no  information  about  the  nhase 
or  amplitude  spectra  of  the  signals,  which  are  statistical  variables,  but  de- 
pends on  quantities  which  are  true  constants  of  the  motion.  For  exaniple,  it 
is  more  certain  that  the  phase  angle  between  vertical  and  radial  components 
is  90  degrees  for  Rayleigh  waves  than  that  the  oower  spectrum  takes  on  a 
given  value  at  20-second  period  for  events  from  Kamchatka. 
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GAIN  VERSUS  INPUT  SNR  FOR  TCA 
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SECTION  V 


CONCLUSIONS 

In  this  report  the  performance  of  three  processors  for  single- 
site  seismic  data  is  presented.  The  processors  are  the  Wiener  filter,  the 
matched  filter,  and  the  three  component  adantive  filter.  Signals  were  buried 
in  noise  at  known  signal-to-noise  ratios,  and  the  outputs  of  the  processors 
were  comoared  to  that  of  a bandpass  filter. 

The  Wiener  filter  used  a synthetic  reference  power  spectrum 
and  had  its  greatest  improvement  in  signal-to-noise  ratio  when  correlation 
terms  between  the  signal  and  noise  were  ignored.  The  gain  of  this  filter  over 
the  bandpass  filter  ranged  from  2 to  6 dB,  depending  on  the  noise  and  signal 
samples  examined, 

A prewhitened  matched  filter  using  the  same  reference  event 
(as  the  Wiener  filter)  displayed  from  6 to  10  dB  improvement  over  the  band- 
pass filter.  By  the  nature  of  the  matched  and  Wiener  filters,  the  matched 
filter  is  expected  to  have  better  gain. 

The  three  component  adaptive  filter  yielded  from  1 0 to  17  dB 
improvement  over  the  bandpass  filter  when  a weighting  function  based  on  the 
probability  distribution  function  of  the  variables  of  the  surface  wave  particle 
motion  was  employed. 

1 

t 

For  all  processors  the  gains  found  when  signals  were  held  con- 
^ stant  were  more  consistent  than  when  noise  samples  were  held  constant.  This 

is  presumably  because  details  of  the  noise  are  less  important  in  determining 
gain  than  are  details  of  the  signal  being  sought. 
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